Density of proper delay times in chaotic and integrable quantum billiards 
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We calculate the density P(t) of the eigenvalues of the Wigner-Smith time delay matrix for two- 
dimensional rectangular and circular billiards with one opening. For long times, the density of these 
so-called "proper delay times" decays algebraically, in contradistinction to chaotic quantum billiards 
for which P(r) exhibits a long-time cut-off. 

PACS numbers 05.45.Mt, 42.25.Bs, 73.23.-b 



In classical mechanics, the length of time that a sin- 
gle particle remains within a certain region of space 
is a uniquely defined quantity. In quantum mechan- 
ics, however, the concept of "time spent in a region" 
is not well defined, as there is no such thing as a Her- 
mitian "time-delay operator" .El For scattering in one di- 
mension, Wigner has constructed a quantum-mechanical 
delay time in terms of the energy-derivative of the phase 
shift acquired upon scattering from the region or potent 
tial of interest.El His concept was generalized by SmithHI 
who introduced a matrix of delay times, 
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where S is the scattering matrix, s is the energy of the 
incident particle, and h is Planck's constant. In addition 
to being relevant for the retardation of a wave packet, 
the Wigner-Smith time-delay matrix Q has been .shown 
to be related to the capacitance,QB thermopowerR and, 
indirectly, to parametric conductance derivatives!!] and 
quantum pumping.oc3 

In this note, we consider the eigenvalues of the Wigner- 
Smith time-delay matrix for scattering from a cavity. The 
eigenvalues-Tn , 1, . . . , N of Q are known as "proper de- 
lay times" n A cavity coupled to the outside world via 
a waveguide with N propagating channels at energy e 
is characterized by N proper delay times. In the semi- 
classical limit of large N, the system is described by the 
density P(t) of delay times. Our aim is to compare P(t) 
for the cases where the classical dynamics of the cavity 
is chaotic or integrable. On the level of classical delay 
times, it is well known that the delay time distributions 
are different for these two cases: For an integrable cavity, 
Pc\ass(i~) oc t~ 7 has algebraic tails for large r, whereas 
Pciass decays exponentially—foii large t if the cavity has 
chaotic classical dynamics.EJtj Simple arguments based 
on the proximity of initial conditions to trapped periodic 
orbits set 7 < 3 for the two dimensional square and cir- 
cular billiards, the two examples of integrable cavities we 
consider here. 

The density -P(t) of proper delay times for a chaotic 
cavity was caLculated by Frahm, Beenakker, and one of 
the authora!3'L3 using random matrix theory.Ej It was 
found that P(r) has a finite support, 



P(t) = I 27T 



T_ < T < T4 

otherwise, 



(2) 



where t± — f (3 ± v°)j t being the average delay time. 
The density of proper delay times for an integrable cavity 
is markedly different, as we now show by consideration of 
the rectangular and circular billiards in two dimensions. 
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FIG. 1. Density P of proper delay times for a rectangu- 
lar billiard (shown in inset). The density P is normalized to 
unity, J P(r)dr = 1, and the delay times r are measured in 
units of their average f . The inset shows the long-time tail of 
the cumulative distribution P(r')dr', together with a (dis- 
placed) linear fit with a slope of -1.70. The delay times were 
taken for an incident wave with energy near the 17,600th level 
of the billiard, at which the contact has N = 41 propagating 
modes, and the density shown was obtained after averaging 
over small size-preserving variations of the shape of the bil- 
liard and the position of the lead. The circles denote the cu- 
mulative density for a single realization near the 3.5 x 10 6 th 
energy level in the billiard with iV = 585 and a fit slope of 
-1.66. 



To find the density of proper delay times for the inte- 
grable billiards, we solve the Schrodinger equation in the 
billiard and the leads separately, andjaatch the wave- 
functions at the billiard-lead interface.113 Repeating this 
process for all possible incident modes in the lead, a sys- 
tem of linear equations is formed from which the scatter- 
ing matrix S, its derivative dS/de, and, hence, the time- 
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FIG. 2. Density of proper delay times for the circle billiard 
(shown in inset) and the corresponding long-time tail of the 
cumulative density J P(r')rfr'. The distribution shown is 
obtained after averaging over small variations of the energy 
around the 2250th energy level of the billiard, for which there 
were 18 conducting modes in the point contact. 



delay matrix Q can be calculated. The billiards and the 
leads are shown to scale in the insets of Figs. |l| and |^. 
To improve our statistics, we have performed an average 
over the position of the lead and small area-preserving 
fluctuations of the billiard aspect ratio for the rectangu- 
lar billiard, and the small fluctuations of the energy e 
for the circular billiard. Plots of the ensemble-averaged 
P(t) for these two billiard are shown in Figs. |l| and ||, 
respectively. We have studied the long-time asymptotics 
through the integrated density f P(r')dT', see Figs, [j] 
and |^, inset. For large N, the averaged density is rep- 
resentative of the density P(t) of a particular billiard, 
as is shown in Fig. |l|, where we compare the tail of the 
averaged P{t) with the density of proper delay times for 
a single realization. 

As can be seen from Figs, [l] and g, the density P{t) 
was found to decay algebraically for large r for both 
the rectangular and circular billiards. In the rectan- 
gular billiard, the exponent of the decay was estimated 
as 7 = 2.6, independent of energy in the inspected en- 
ergy range 10 3 <e<2xl0 4 (energy is measured in 
units such that the level spacing in the closed billiard 
is one), giving a range in the number of propagating 
modes in the contact ranging from N = 9 to N = 41. 
In the circular billiard we inspected energies in the range 
10 3 < e < 3 x 10 3 (corresponding to N ranging from 
14 to 20), and found that P decayed with an algebraic 
tail with exponent 7 = 3.05, again independent of en- 
ergy within our accuracy. (The result that 7 is larger 
in the circular case could be attributed to rounding er- 
rors in calculating Bessel functions in the large N regime 
of interest, introducing small random fluctuations that 
mimic a small randomness in the potential, thereby sup- 
pressing long dwell times as in the chaotic regime. The 
jaggedness of the small r distribution in the circular case 



is due to the restricted avenues of averaging available as 
compared to the rectangular case.) The classical distri- 
bution_pf delay times decays oc t~ 7 with 7 rs 3 in both 
cases.EZl Hence, we conclude that the density of proper 
delay times decays algebraically for the integrable bil- 
liard we studied, with a power that is close, but (for our 
observations) not precisely equal to the classical power 
7 = 3. Our results majz-be compared to one dimensional 
quasi-periodic systemslia which also exhibit an algebraic 
decay of the delay time distribution. 

We wish to point out that the difference in the den- 
sity of delay times for chaotic and integrable cavities 
reported here is something remarkable. Although it 
is known that many quantum properties are different 
for cavities with chaotic and integrable classical dy- 
namics, these differences usually pertain to the statis- 
tical fluctuations, described by correlation functions, or 
show up in small quantum Jjitcirference corrections to a 
(semi)classical background ,EJEfl but not in the (ensem- 
ble averaged) densities themselves. The only exception 
known to the authors is the density of states p(e) in a so- 
called "Andreev quantum billiard" , an "electron billiard" 
that is connected to a superconducting point contact. In 
this case, electrons that exit the cavity and impinge on 
the superconductor interface are reflected as holes and 
vice versa. As a result of thia-special reflection process, 
known as Andreev refiection,E2l p{e) is singular around 
the Fermi energy e = 0: The density of states has a gap 
for a cavity with chaotic classical dynamics, while foiian 
integrable cavity, p{e) decays algebraically as e —> O.EiJ In 
fact, this feature can be connected to the difference for 
the tail of the density of proper delay times P(r) reported 
here via the following heuristic argument: As shown in 
Ref. [2l], the Andreev quantum billiard has an eigenstate 
at energy e precisely if the matrix product S(e)S'(—e) 
has an eigenvalue —1. Expanding S around e — gives 

S(e)Sfl(-e) « e 2leQ ' h . (3) 

With this approximation, the condition for eigenstates in 
a cavity coupled to a superconductor simplifies tcE3 

e 2ier n /h = _ lj ( 4 ) 

where r n is a proper delay time (eigenvalue of Q). 
Eq. ([I])), viewed as a constraint on the product er„, con- 
nects p(e) at small e to P(t) at large r. The fact that 
the density of states is gapped for the chaotic Andreev 
quantum billiard can then be understood as following 
from the absence of a large-time tail of P(t), cf. Eq. (0), 
whereas the algebraic vanishing of p(s) near zero energy 
for the integrable Andreev quantum billiard is seen to be 
related to the. algebraic tail of the density of delay times 
in that caseEj We note, however, that the approximation 
Eq. (^|) cannot be used for a quantitative estimate of the 
density of states, as it becomes unreliable for energies 
er ~ h, which is precisely where the first eigenstates are 
expected to appear. 
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While we are not aware of a method to directly mea- 
sure the density of proper delay times for an electronic 
system, a direct measurement of P(r) would be possi- 
ble using the scattering of electromaffiaetic waves (mi- 
crowave radiation) from a metal cavity.E-3c3 With a suit- 
able choice of-, basis, S and Q may be simultaneously 
diagonalized.tij If the incoming electromagnetic radiation 
is in a plane wave state corresponding to one of these 
basis vectors and slowly modulated in intensity, the ac 
modulation of the outgoing waves will be delayed by the 
propejndelay time r„ corresponding to the incoming wave 
modefl In this way, the qualitative difference between 
chaotic and integrable cavities should be readily acces- 
sible from the tail of the density of proper delay times, 
without further ensemble averaging. 
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dation, and by the Natural Sciences and Engineering Re- 
search Council of Canada. 
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